Phase transition of interacting disordered bosons in one dimension 
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Interacting bosons generically form a superfiuid state. In the presence of disorder it can get 
converted into a compressible Bose glass state. Here we study such transition in one dimension 
at moderate interaction using bosonization and renormalization group techniques. We derive the 
two-loop scaling equations and discuss the phase diagram. We find that the correlation functions 
at the transition are characterized by universal exponents in a finite region around the fixed point. 
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Disordered quantum systems show a remarkable rich- 
ness of phenomena. One of the most spectacular mani- 
festations of the interplay between quantum effects and 
disorder is provided by Anderson localization, where the 
interference between waves can lead to exponentially lo- 
calized wave functions. For noninteracting particles, An- 
derson localization is by now well understood and has 
been spectacularly tested in various experiments ranging 
from microwaves to cold atoms [l[ . The situation is con- 
siderably more complex when interactions are present, 
and a full understanding of disordered interacting quan- 
tum systems is largely still missing. Among such systems, 
particularly interesting are disordered bosons. Indeed in 
the absence of disorder, interacting bosons are expected 
to undergo a Bose-Einstein condensation and to become 
superfiuid (SF), and hence naively should be able to resist 
disorder. One can thus expect the competition between 
disorder and interactions to be particularly relevant in 
that case. In addition since all noninteracting bosons 
collapse into a disorder created lowest energy state, in- 
teractions must be included from the start to understand 
this competition. 

Indeed it was shown by a renormalization group (RG) 
analysis Q that disordered interacting bosons would un- 
dergo a phase transition between a SF and a localized 
phase. This transition and this phase, nicknamed Bose 
glass (BG), was also shown [3j by scaling arguments to 
exist in higher dimensions. In one dimension the SF- 
BG transition is in the Berezinskii-Kosterlitz-Thouless 
(BKT) universality class, and universal exponents at the 
transition exist for the various correlation functions, in 
particular the single particle one. One interesting possi- 
bility, which was argued in Q was the possibility of the 
existence of two localized phases, the BG corresponding 
to a strong interaction, strong disorder fixed point, while 
another would correspond to a weak interaction, strong 
disorder fixed point. These predictions on the phase di- 
agram and the flow were found to be in good agreement 
with numerical studies of this problem jj HJ . 

Recently this problem has regained a considerable in- 
terest thanks to its remarkable realization in cold atomic 
gases [6-8]. On the theoretical side, new studies focused 



on the low interaction, strong disorder case 0-[T3|. In 
particular using a real-space renormalization group study 
of a related disordered Josephson junction array model a 
BKT-like transition was found again, but now the Lut- 
tinger parameter takes disorder dependent value at the 
transition |13j |. These results are thus in contradiction 
with the universal exponent found by the RG analysis Q 
for intermediate interactions, which would then suggest 
the existence of two distinct phases. It is thus important 
to ascertain that the universality of the exponent found 
in [2| is not an artefact of the use of the lowest order RG 
but survives if higher order are taken into account. This 
is important in view of what can happen for other disor- 
dered systems such as the Cardy-Ostlund model [3] for 
which some nonuniversal terms are generated at next or- 
der. In addition pushing the RG flow to the next order is 
also useful in view of the comparison with e.g. numerical 
studies of the phase diagram. 

In this letter we perform an analysis of the SF-BG 
transition using a Luttinger liquid (LL) description of the 
interacting bosons and find an RG flow to second order 
in disorder. We show, both from the RG flow and from 
symmetry arguments that the exponents remain univer- 
sal at the transition whenever the LL description can be 
used. This implies that the exponents will stay univer- 
sal on at least a part of the phase diagram. If there are 
indeed varying exponents at weak interactions a transi- 
tion between two localized phases must exist. We also 
clarify the respective influences of the forward and back- 
ward scattering on the impurities and compute the single 
particle and density correlation functions. 

We describe a quantum interacting Bose fluid in one 
dimension. Its low energy properties can be described by 
the Tomonaga Luttinger liquid model [IH which corre- 
sponds to the Hamiltonian 



H = ±-Jdx [vK[d x e{ x )] 2 + -^[a^(.*)] 2 } , 



(i) 



where v denotes the sound velocity and K is a dimen- 
sionless Luttinger parameter. The fields and ip satisfy 
[6(x),dytp(y)] = iirS(x-y). The density of particles reads 
p(x) = po + 2pid x ip(x) + 2/?2 cos[2,tp{x) — 2irpQx], where 
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pi = — l/2ir, p is the average density, and P2 is a nonuni- 
versal constant. For a local repulsion q between bosons, 
one has < 1/K < 1 for < g < oo. [HI 

In addition we put the system in a disordered potential: 



H d = J Ax {2p lV (x)d x <p + p 2 [C(x)e i2 * + h.c] } , 



(2) 



where we distinguish the forward 77 and backward £ 
scattering part of disorder with Fourier components 
around momenta and ±27rpo, respectively 15]. We 
assume Gaussian disorder with correlations rj(x)rj(x') = 
h 2 D f S{x - x'), £(x)£*(x') = h 2 D b S{x - a/), while the 
other correlations vanish. By ... we denote the dis- 
order average. The total Hamiltonian of the system is 
H = H + H d . 

We average over disorder using the replica trick. The 
replicated Euclidean action is S = Sq + Sf + Sb- The 
quadratic part reads 



So 
h 



2irK 



1 



(d x (p a ) + -^(d T ip a ) +m (<p a ) , (3) 



where ip a (x, r) is a set of bosonic fields. Here a — 1 ... n 
and the limit n — > has to be taken in the end. We intro- 
duced the shorthand notation J XT . . . = ^ Q / dxdr . . ., 
where Greek letters denote replica indices. As an infrared 
cutoff we introduced a small mass m which will be sent 
to zero at the end of calculations. The disorder part of 
the replicated action reads 

Sj_ 

h 

Sb 

Jl 



: -2p 2 D f I [d xVa {x,r)][d xV p{x.r% (4) 
-p%D b I / cos[2cp a (x,T) - 2ipp(x,T')]. (5) 



Introducing the Fourier transform one can diagonalize 
the harmonic terms ([3]) and ((4]) in S. The correlation 
function G a p(x,T) — {ip a (x, r)tpp(0, 0)) is then written 
17] as G a p — G5 a p + Gf, with 



G(x,t) = —Ko{m\/ x 2 + v 2 t 2 + a 2 ), 



Gf(x) 



TT 2 K 2 pjDf 



■1*1(1 



m\x\ 



(6) 
(7) 



Here the parameter a has been introduced as the ultra- 
violet cutoff, and (...) is the average with respect to the 
harmonic part (So + Sf)/h. Kq denotes the modified 
Bessel function of the second kind [18j . 

To treat the backward scattering part of the disor- 
der, we use a field theoretic approach 1^-21| to ob- 
tain the effective action V of the model. For an ac- 
tion S(cp), it is defined as T(ip) = Jtp — W(J), where 
W (J) is the generator of connected correlations given by 
W(J) = In/ Vipe- S ^^ h+Jv . Using J(x) = j^y, we ob- 
tain T{(p) = - In / V X t 



-S{ V +x)/h+fdx X (x) 



til") . This 



general equation can be solved perturbatively, order by 
order with respect to the small parameter Db filf . The 
perturbative expansion of the effective action contains all 
the information about critical properties of our model (at 
the first two orders in Db, in our case). In order to derive 
the scaling equations one should find the divergent terms 
in r in the limit a — > 0. After expanding the operators 
one finally obtains the effective action [22j 



T=E/ dxdr^[{d xVa ) 



m 2 (ip a ) 2 



1 



2ttKv 



2.3§ ai + 0{38 2 ) 



(d T ip a ) 



[2p 2 D f + 2^ 2 b 2 ] 



^2 / dxdTdT'[d x tp a (x,T)}{d x ipi3{x,T r )} - [& + 2b 1 .@ 2 ] ^ / dxdrdr' cos[2(^ Q (a;, r) - 2(fp(x,T% (8) 

a/3 J a f) •* 



where we have introduced 3B — p\Db<& 4G (°'°). The co- 
efficients ai,6i, and b 2 read a\ — / drr 2 [e 4G (°' T ) — l] , 
b 2 = J dxdrdr'x 2 ^ 40 ^^ - l][e iG ^ - 1], and 
61 = / dxdrdr' f 2 (x,T,r + t'), where f 2 (x,T,r') = 

e 4G(i,r)-4G(i,T')+4G(0,T-T') _ & iG(x,T) _ e -4G(x,r') _ 
e 4G(0,r-r') + g + 4 [ e 4G(0,T-r') _ l][Q(x, t') - G(x, t)]. 



We can now calculate the unknown parameters in (jHJ . 
Introducing the small parameter 5 = K — 3/2 which mea- 
sures the distance from the critical point (see below) and 



using © we get 



-2A + SX 2 - 4(ln2 - 1)SX + 0{5 2 ) + 2c x 



2 (emu) 3 



(9) 



where A = \nc 2 m 2 a 2 , c 



= e 7B /2, 7s is the Euler con- 
stant, and Cj is a constant. The next two terms are 

7T [9A 2 + 2(51 - 54 In 2) A + 0{5) + 4c 2 ] 

°1 = o..9/„-^ ) ( 10 ) 



8v 2 (cm) 3 



b 2 = 



;-.V+0(6)} 



v 2 (cm) 6 a 



(11) 
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The effective action (8) has divergencies when ma — > 
contained in A. In order to remove them we introduce a 
set of renormalized coupling constants (denoted by the 
subscript R ) by 9 = Z h 9 R , v = Zv R , 8 = Z(3/2 + S R ) - 
3/2, K = Z(3/2 + 5 R ),m R = m where 

Z b =1 - - [(39 - 54 In 2 + 9ci) @ R + 2S R ] A 

+ ~(9® R + 25 2 R )\ 2 ~3> R (6 Cl + c 2 ), (12) 



Z =1 - [3£r + (6 In 2 - 4) 9 R S R ) A + -^feA 2 
+ 3ci^ fl + 2ci9 R S R . 



(13) 



We introduced the dimensionless disorder strength 3 = 
irp 2 a 3 Db/v 2 . The effective action expressed in terms of 
renormalized quantities does not contain divergencies. 
This leads to the RG equations of the coupling constants, 
obtained by requiring that derivatives of the initial cou- 
pling constants with respect to m R nullify. They read 



d£ 
dSR 
di 

^D fR 



= -2® R 5 R + A% + 0{9%6 R ), 
= -9@ R + BS> R 6 R + G(@ 2 R ), 
ii. — l-rM=0, 



VR 
Kr 



(14) 
(15) 
(16) 



where l~— \nrriR. The constants A = 54 In 2 — 39 — 9ci 
and B = 6 — 18 In 2 + 9ci depend on the details of the 
choice of regularization in ([3j 23]. However, their sum 



A + B = 361n2 - 33 



(17) 



is a universal number for the model. Equations (|14|) -(fl7 |) 
are the main results of this letter. 

At the first order in the backward scattering (i.e., 
putting A = B = 0) the above equations are identical 
to those first derived in Giamarchi and Schulz @ (see 
also [HI). One could notice the absence of renormaliza- 
tion to v/K which an exact result in all orders due to 
the statistical symmetry <p(x, r ) — > tp(x, r) + w(x) of the 
disordered part of the action [24j, [25j. The second or- 
der terms are novel and lead to additional insight in the 
physics of the SF-BG transition for this model. At the 
fist order in Db no off-diagonal (in replica space) renor- 
malization of the quadratic action could be generated, 
for trivial reasons. The second order is the lowest order 
at which such terms might in principle appear. From 
([8]) we see that no such terms are generated, except for 
terms describing the forward scattering process, which 
have no consequences for the localization properties of 
the system. In addition, the forward scattering terms 
contain no divergence and hence there is no renormal- 
ization of Df. The quadratic part of the effective action 
thus remains, even to the order D 2 essentially diagonal 
in replica. The transition between the SF and the BG 



phase can thus be fully characterized by the two param- 
eters K and Df, and not, like in other models such as 
the classical Cardy-Ostlund model [3], by a full set of 
variables (corresponding to the off-diagonal terms) that 
enter into the correlation functions and can affect the ex- 
ponents of correlation functions at the transition. This 
important result, proven here directly from the RG flow, 
is a consequence of the time independence of the disorder 
[26j ]. Indeed for two independent replicas a ^ j3 before 
averaging over disorder one obtains 

(<f a (x,T)(p (x',T')) =(<p a (x,T))(<pp(x' ,t')) 

=(ip a {x,Q)){ip p ^,0)) (18) 

since the disorder does not depend on r. The correla- 
tion (|18|) is thus time independent. Thus terms such 
as J dxdx'dT[d T ip a (x, t)] [d T (pp(x', r)] cannot appear in 
the effective action, at any order, since they would lead 
to time dependence for (|18[) - This has important con- 
sequences for the physical properties at the separatrix 
between the phase for which the disorder is irrelevant 
(Db —> 0) and the phase for which the disorder is relevant 
(the BG phase). The absence of such off-diagonal replica 
terms thus leads to a universal value of the parameter 
K, namely K c = 3/2, and correlation functions will thus 
decay with a universal exponent at the transition. Our 
analysis thus confirms that the SF-BG transition around 
the value K = 3/2, i.e. for intermediate interactions, has 
a generic universal exponent in a finite region around 
that point. This puts stringent constraints on the phase 
diagram, as schematized in Fig. Q] Indeed this result 
that uses the bosonization representation of the disor- 
dered Bose gas is guaranteed to work when the disorder 
is weaker than the chemical potential. This means that 
for a good part of the SF-BG boundary the exponent 
will remain constant. What happens for larger disorder 
or weaker interactions is still an open question. One pos- 
sibility is of course that the exponent remains universal 
along the whole line. In that case, there is most likely 
a unique BG phase. The only singular line of the phase 
diagram would be in that case, the non-interacting line 
U = for which the bosons are localized by rare events 
of the random potential in a finite region of space. How- 
ever if for weak interactions one can obtain non- universal 
exponents at the transition SF-BG as discussed in [l3j ]. 
it is then impossible to smoothly connect the two sepa- 
ratrices between the SF and the BG. It would thus imply 
that there is on the separatrix a critical point above which 
the exponents would be universal and below which they 
would vary with parameters. This would have for conse- 
quence to have two different localized BG phases, as was 
argued as one of the possible scenarios in [2j . A possible 
order parameter distinguishing these two localized phases 
might be the moments of SF stiffness distributions. An 
interesting question is whether a similar mechanism oc- 
curs inside the superfluid as well (see e.g. 27[). 
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The solution of the flow equations ([14]) and (|15l) is 

5% - 9$> R - A5 r $>r + 2(A + B)8%/27 = C, (19) 

where C is an arbitrary constant, and C < (C > 0) 
marks the insulating (SF) phase. For C = the sys- 
tem is at the critical line, where &r = (Sr/3) 2 + (2B — 
A)S 3 R /2A3 + 0{8 R ). Then we obtain the solution of (H3 
at large scales 



B) hvi 



(20) 



FIG. 1. (color online) Possible phase diagram of a disordered 
ID Bose gas as a function of the boson repulsion U and dis- 
order S>. The green dash-dotteed line schematically indicates 
the boundary below which a bosonization description of such 
a problem is guaranteed (namely the disorder is smaller than 
the chemical potential). In this region the superfluid-Bose 
glass transition, denoted by the solid blue line, is described by 
universal exponents. The question of the nature and critical 
behavior of the transition in the regime for which bosoniza- 
tion cannot be directly applied is yet open (see text). One 
possibility (not shown on the figure) is that the exponent 
remains universal along the whole SF-BG line, leading to a 
single BG phase. The only singular line is then the noninter- 
acting bosons (U = 0) which are localized by rare events of the 
disorder and indicated by the blue box. The other possibility 
(shown in the figure), is that if at small interaction, there are 
nonuniversal exponents along the dashed blue line, then this 
forces the existence of a critical point on the SF-BG boundary 
between the nonuniversal and universal regime. This would 
imply the existence of two distinct BG phases. 



In general the solution of RG equations ([19]) is nonuni- 
versal. However along the critical line 5r(£) is universal 
since the invariant of the model (]17l) emerges, which in 
turn determines universal form of the correlation func- 
tions at the transition, as shown below. Since 2>r ~ 8 2 R 
close to the transition the neglected term 0{3§ 2 ) in (|5J) 
would produce 0(@ R ) in ([15]) which is beyond the two 
loop order we consider. The correlation length close to 
the critical line from the insulating side (C < 0) takes 
the form £ oc exp(n/ y/9® R - 6 R )[l + 0(8 R )}. Therefore 
higher order corrections do not affect £ in an essential 
way and the transition is of BKT type [2j . 

The two loop RG also allows us to compute the sin- 
gle particle correlation function F(x) — (x)^ (Q)) h 
at the transition with improved accuracy. It can 
be obtained from the Euclidean action S as F(x) = 
p ^ e i[8(x,o)-e(o,o)]^ s _ Performing the perturbation theory 
with respect to ([5]) and taking into account the flow of 
parameters one obtains for \x\ a 



The SF-BG transition can be seen directly from the 
supcrfluid correlation function, which will go from a di- 
vergent power law behavior in the SF phase, to an ex- 
ponentially decreasing one in the localized phase. Vary- 
ing the strength of the disorder and of the interactions 
allows to probe the universality of the exponent at the 
transition. Such behavior can be probed in cold atomic 
systems or in magnetic insulators 28-31]. Cold atoms 



offer the advantage of the control over disorder and in- 
teractions. However optical lattice systems suffer from 
the inhomogeneity due to the confining potential, which 
complicates the analysis of the exponent. Atom chips 
realizations are up to now limited to relatively small in- 
teractions (A" ~ 40). Magnetic insulators are very ho- 
mogeneous and allow precise control and measurement 
of the boson density and compressibility. Controlling the 
disorder is, however, more challenging. No doubt that 
for both types of systems these difficulties will be over- 
come in the future, allowing reliable answers to the above 
points. 

Let us now use our improved RG flow to investigate 
more quantitatively the behavior close to the transition. 



F(x) = po 



1/3 



111 - 



2/9 



1 - 



2(A + B) In In 



248 



hi 



kl 



(21) 



Note that forward scattering term does not appear in 
F(x). In addition to the lowest order result F(x) oc 
| | 1 / 3 Q j W e find a logarithmic correction ln 2 ^ 9 (|x|/a). 
Finally the role of the universal quantity (]17p is found in 
subleading corrections to F(x) in (|2"T]) . 

The 27rpo part of the density-density correlation func- 
tion F{x,t) = (e l2 <^. T ) e - l2l ^ < ))s can be obtained in a 
similar way. On the separatrix it reads 



F{x, t) oc exp 



^ 2 ^p\D f \x\ 



2{A + B) lnlnf 



27 



In- 



(22) 



where R = yjx 2 + v r t 2 » a. The first term of (|22J is 
due to the forward scattering term, and we see that is 
comes with DfK 2 /v 2 that is not renormalized. 

In conclusion, we have computed to two loop order the 
RG equations for interacting disordered bosons. We have 



shown, both from the RG and from general arguments 
that in the regime described by bosonization (intermedi- 
ate interactions) the separatrix between the BG and SF 
is characterized by a universal exponent. We computed 
the single particle and density-density correlations at the 
transition. Our calculation uses bosonization and there- 
fore also applies to disordered fermions. Detailed aspects 
of that and the details of calculation will be published 
elsewhere. 
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The defining equation for the effective action from the main text reads 

T(ip) = - In J Vxe~ S{v+x),h+S Axx{x) ^) . (23) 

This general equation can be solved perturbatively, order by order with respect to the small parameter [l| . Up to 
an additive constant, we obtain r = (So + Sf)/h + T 1 + T 2 + O(Df), with Ti = (SbQp + x))x/^ anci 

ra _ _ + 1 r; + 1 zf^r'C^ -,,r- r'^^ry (24) 

where (• • • ) x denotes an average with respect to the field x- From the technical point of view, our approach is very 
similar to the standard Wilsonian approach when one performs thinning of high-momentum degrees of freedom. In 
the present case the difference is that the fields \ that are integrated out live in the whole momentum space and in a 
cumulant expansion calculates the effective Hamiltonian. Averages that appear during such procedure of integrating 
out high-momentum degrees of freedom are with respect to the quadratic Hamiltonian and one uses Wick's theorem, 
see for example Ref. Q for details. 
At the lowest order we easily obtain 

Y x =-m s ^ l [ dxdrdT'| e 4G (°< T ~ T ') - 1 8 a p + l} cos[2ip a (x, r) - 2<p fi (x, t% (25) 

where S3 = p?,Dbe~ 4G ( ' r '\ One should notice that only a part of the full correlation function that is diagonal in 
replica indices, (6), enters the result. This is due to the fact that the off-diagonal part (7) of the propagator does not 
depend on imaginary time. 

The second order term is more complicated and could be written as a sum T 2 = Sj=i ^2 s where j denotes the 
number of replica sums in the expressions: 

if' =-\sS 2 ( f 1 (x-x 1 ,T,T\n,T{)cos[2<p a (x,T)-2cp a (x,T , ) + 2ip a ( Xl> n)-2 l p a ( Xl> T{)], (26) 



2 



C 1 T 1 T 1 



T { £ ] =-£g 2S p I \ f 2 (x- xi,t -t{,t -Ti)cos[2(p a (x,r) +2(p a (xi,Ti) -2(p a (xi,r{) -2(pp(x,T')} 

S = ±l J a ff k 

+ \hi x ~Xi,T - n,T - t[,s) cos{2tp Q (x,T) - 2sip a (x 1 ,T 1 ) - 2ipp(x,T') + 2sipp (xi,t[)] I , (27) 
r^ 3) = - SS 2 S~] / fi( x - Xi,t - ti, s)cos[2(p a (x,T') - 2ip fj (x, r) + 2sipp(x 1 , n) - 2s^(xx,t[)]. (28) 



S=±l a p-t 

The / functions in the previous equations are / 2 (x,t,t') = e ^G(x,T)^4G(x,r')+4G(o,T-r') _ e 4G(x,r) _ e -4G(x,r') _ 

e 4G(0,r-r') + 2 + 4 [ e 4G(0,T-r') _ \][G{x,t') - G( X ,t)], / 3 (x,T,t',s) = [e isG ^ - 1] [e isG ^ - 1], f 4 (x, T, s) = 

e 4sG(x,r) _ 4 s c( X T ) _ i The expression (121)1) leads to corrections of K of the order 0(D 2 ), that is a higher order in 
the expansion than the one considered in Eqs. (15) and (19) of the main text . 

The obtained perturbative expansion of the effective action contains all the information about critical properties of 
our model at the first two orders in D\,. In order to derive the scaling equations one should find the divergent terms 
in r in the limit a — > 0. After expanding the operators one obtains the effective action 

T=J2 J dxdT {^ [(d x <Pa) 2 +m 2 ( Va ) 2 ] + ^-L- + 233 ai +0(3§ 2 ) (c^ Q ) 2 j - [2p\D } + 2SS 2 b 2 \ 

x / d%dTdT'[d x <p a (x,T)][d x <pp(x,T')} - [38 + 2b 1 S8 2 ] ^ / dxdrdr' cos[2<^ Q (a;, r) -2^(x,r )]. (29) 

a/3 a/3 J 
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The coefficients ax, bx, and b 2 arise from the expansion of operators in ([25]) and (|2"7|) and read ax — J drr 2 [e 4G( -°' T ' — l] , 
bi = f dxdrdr' f 2 (x,r,r + r'), and b 2 — J dxdrdr' x 2 f 3 (x , r, r', 1). The terms from (j27|) produce corrections to 
operators (8), contained in b\ and b 2 . One can notice that due to the parity / 3 (x, r, r',1) = / 3 (x, — r, — t', 1) the 
analogous term to (4) but with time derivatives is absent in (8). This differs from what happens in the the related 
classical 2D Cardy-Ostlund model [3( where such terms are generated. The three-replica part (|2"5)) contains two 
contributions that turn out to be unimportant. Namely, for s = 1 after a gradient expansion a free sum over /3 
index delivers a factor n, that vanishes in the replica limit. For s = — 1, one generates an operator of the form 
oc cos{2(p a + 2<p 7 — 4^3^) which has a nondivergent prefactor and in addition is irrelevant close to the critical point. 
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